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ABSTRACT 

We present a new approach to the analysis of time symmetry in light curves, such as those in the 
x-ray at the center of the Scorpius X-1 occultation debate. Our method uses a new parameterization 
for such events (the bilogistic event profile) and provides a clear, physically relevant characterization 
of each event's key features. We also demonstrate a Markov Chain Monte Carlo algorithm to carry out 
this analysis, including a novel independence chain configuration for the estimation of each event's 
location in the ligh t curve. These tools are applied to the Scorpius X-1 light curves presented in 
I Chang et al.l ()2007f ). providing additional evidence based on the time series that the events detected 
thus far are most likely not occultations by TNOs. 

Subject headings: Kuiper Belt - methods: data analysis - methods: numerical - methods: statistical 
- stars: individual (Scorpius X-1) 
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1. INTRODUCTION 

As new surveys acc umulate vast databa s es of astro- 
nomical light curves (|Hodapp et al.l 12004 iStarr et all 
[2002) . the identification and analysis of events in these 
time series is becoming a n increasingly im portant task 
(e.g. iPreston et all (|2009f) : iStewartI (|2009D ). One of the 
topics of such analyses is examining the symmetry of 
such events. Traditionally, symmetry has been analyzed 
by examining the skewness of each event, so metimes after 
stacking the light curves (for examples, see iJones et al.l 
(|2006l l2008n ). However, this approach is not satisfac- 
tory for x-ray light curves due to the small number of 
counts involved, particularly during suspected occulta- 
tion events. Additionally, the nonparametric approach 
of examining skewness has significantly lower power than 
a parametric analysis with a well parameterized model, 
especially considering the large r nean squared errors of 
common skewness estimators (see IJoanes fc GillfTOOl ). 
This means that, for a given false-positive rate, a well- 
designed parametric analysis will have a lower false- 
negative rate than a nonparametric approach. 

The debate over the nature of the observed Scor- 
pius X-1 events has raged in the literature since 2006 . 
It began with the pubhcation of iChang et al.l (|2006D . 
which claimed the detection of at least 58 occultations 
of Scorpius X-1 by small (diameter < 100m) TNOs. 
The size distribution of these small TNOs bears impor- 
tant clues of th e dynamical evolu ti on of the early So- 
lar Sy stem (e.g. IPan fc Sanl (|2005l ): iKenvon fc BromlevI 
(|2004f )). yet the vast majority of TNO population is still 
way beyond the limit of direct observation by even the 
largest ground-based telescopes. Occultation of back- 
ground stars by TNOs is the only present observational 
method able to reach the smaller obje cts (Bailcv 1976; 
ILehner et ani2008t IRoques et all 120061: iBickerton et al.i 
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120081: IBianco et alll2009f). Th i s find ing was challenged 
later that vear by IJones et al.l (j2006f ). The latter group 
claimed that the observed events were unlikely to be oc- 
cultations by TNOs because the stacked profile for these 
events was asymmetric. However, this analysis was not 
entirely satisfactory from a statistical perspective; it was 
performed on heavily transformed versions of the origi- 
nally data, did not take into account the discrete nature 
of the observations, and provided no measures of uncer- 
tainty. 

The response of Chang et al. to this analysis was 
IChang et al.l (|2007f ). wherein they claimed to have iden- 
tified 12 events in the Scorpius X-1 data that were un- 
likely to have originated from instrumental e ffects. How- 
ever, this claim was challenged again by IJones etldl 
(|2008f ). Jones et al. presented evidence that the ob- 
served events were due to charged particle events in the 
detectors. Their argument relied again, in part, on the 
claim that the observed event were not temporally sym- 
metric. Unfortunately, their analysis had not proceeded 
beyond a simple s tacking-based approach (analogous to 
IJones et al.l (|2006f )l. so the validity of these conclusions 
remains uncertain^ 

Most recently, ILiu et all (|2008l ) presented a new anal- 
ysis of 72 -ks of data of Sco-Xl taken in the year 2007. 
ILiu et aTl concluded that no significant dips which might 
be real occultation by 60 — 100 m TNOs were observed. 

In this paper we introduce a new parameterization for 
unimodal events in light curves, called the bilogistic event 
profile. This parameterization can describe a wide range 
of possible events and naturally yields estimates of phys- 
ically interesting quantities such as the full width at half 
maximum (FWHM) of the event. Additionally, our pa- 
rameterization can be applied in both the Poisson and 
Gaussian regimes. 

Additionally, we demonstrate a computational ap- 
proach to inference with the bilogistic event profile using 
Markov chain Monte Carlo (MCMC) methods. This in- 
cludes the introduction of a novel approach to the con- 
struction of a proposal distribution for the center of the 
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event, improving the efRciency of our simulations.** 

Finally, we apply our parameterization and computa- 
tional approach to occultation eve nts iden ti fied i n the 
Scorpius X-1 RXT E data in Chan g et alJ (|2006f ) and 
IChang et al.1 (j2007l ). Our analysis lends further support 
to the argument that the events ident ified thus far are in - 
stru mental in or i gin, a s suggested in iJones et"aLl ()2006D 
and lJones et~an (|2008[ ). 

In § 2 we describe the details of our data model for 
the Scorpius X-1 . We present the details of our MCMC 
approach in |§ 3[ The results of apply ing our method to 
simulated data are discussed in |§ 4[ Our primary sci- 
entific results on th e Sco rpius X-1 light curves comprise 
51 We conclude in IS 6 
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2. MODEL 

2.1. Context & Specification 

We have a time series of counts {yt},t £ {0, • • • , T—1}. 
We assume that each count is distributed as an indepen- 
dent Poisson random variable, conditional on the source 
intensity at time t: 

Dt ^ Poisson{Xt). 

We also believe that an event of some kind has occurred 
in this time series. For the sake of this discussion, we 
assume that this event is characterized by a dimming of 
the source (e.g., a transit). However, this framework can 
be extended trivially to events character by increased in- 
tensity (e.g., supernovae) or the high-count regime (with 
a Gaussian distribution replacing the Poisson). The key 
constraint is that we know a priori the sign of the ex- 
pected deviation from the baseline intensity. We charac- 
terize the intensity At as 



\t = 



ag{t:T, 9 ) 



where limj^oo 5(i; t, ) = linit^_oo 5(i; 9 ) = Q and 

supR5(t;r, 9 ) = 5(r;T, 9)^1. We will call g{t]T, 9 ) 
the event profile, as it characterizes the pattern of in- 
tensity changes throughout the event. The parameter r 
is the time index about which the event is hypothesized 

to be time-symmetric, and ^ is a vector of parameters 
characterizing the event, c sets the baseline intensity 
of the source, and a sets the peak size of the deviation 
from the baseline intensity during the event. All of these 
parameters must be estimated from the data. 
Putting all of these points together with the log- 
likelihood for the model, we obtain: 

— + 

£(c, a, T, 6* ; y*) = ^ [yt In (Ai) - Af] -I- const (1) 
t=o 

Following Bayesian approach we can express the poste- 
rior probability as a function of the likelihood and the 
prior probabilities as: 

^ In Mctropolis-Hastings-type MCMC methods, each step of the 
Markov chain consists of two stages. In the first, a new value for 
the variable is generated from a proposal distribution with a known 
density. Then the proposed value is accepted or rejected as the next 
value of the chain with a certain probability (calculated using the 
target and proposal densities). By constructing a proposal distri- 
bution that is similar to the target distribution we are attempting 
to draw from with our MCMC method, we are able to improve the 
efficiency of our algorithm. 



Having obtained the posterior of our parameters (up to 
a normalizing constant), we have all of the information 
needed to make inferences about the parameters of in- 
terest. We can characterize our inferences about the 
locations of parameters via their posterior means and 
medians. Our uncertainty about these parameters can 
be characterized by 68% posterior intervals or posterior 
standard deviations. However, it still remains to specify 

a functional form for the event profile g{t;T, 9). 

2.2. Bilogistic event profile 

The choice of event profile g{t] t, 9 ) is quite impor- 
tant to this analysis. One must achieve the proper bal- 
ance between parsimony in parameterization and giving 
the model enough flexibility to fit a wide range of event 
shapes. To that end, we introduce the bilogistic event 
profile: 



'l{t;T,hi,h2,ki,k2) = 

ht 
kt 



1 + e >'t 



1 + e 

hi t < T 

h2 t>T 

ki t < T 

k2 t>T 



(3) 

(4) 
(5) 



This event profile is characterized by four parameters 
(in addition to the location parameter r). The first two 
parameters [hi and /i2) characterize the (approximate) 
half-maxima of the event's amplitude, represented as de- 
viations from the event's center; we will refer to them 
as the "half-life" parameters. This holds only approxi- 
mately because of the correction for continuity in the nu- 

merator (the e '=* term), which ensures that the event 
profile is continuous at time r. Thus, {hi — ^12) char- 
acterizes the departure of our event from symmetry in 
terms of the half- maxima, and {hi -\- /12) characterize the 
(approximate) full width at half maximum for our event. 

The second set of parameters (fci and fc2) characterize 
the rate at which the intensity changes during our event; 
we will refer to them as the "curvature" parameters. The 
first (fci) characterizes how rapidly our source's intensity 
diminishes at the beginning of our event, and the sec- 
ond characterizes how rapidly the intensity returns to 
the baseline at the end of our event. If fci << hi and 
fc2 << /i2, as is often the case, the curvature parameters 
can been viewed as pivoting our event profile through the 
half- intensity points set by hi and /12, with higher values 
corresponding to sharper changes in intensity. These fea- 
tures of the bilogistic event profile can be seen in Figures 
[1] and [21 where we fix one set of parameters {hi and /i2 
or ki and ^2) and vary the other. 

To perform inference using this model, we first consider 
it in a Bayesian light, assuming flat priors on all param- 
eters except r. We place one restriction on r a priori. 
We assume that the event is centered within the inte- 
rior of our light curve. In this context, we take interior 
to mean far enough from the boundaries that the event 
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is contained within our light curve. This assumption is 
reasonable for our analysis because the light curves we 
are working with have been extracted and preprocessed 
to contain one complete event each. Thus, we impose a 
prior probability of zero on values of r within m millisec- 
onds of the edge of the light curve for our analysis; for 
the analysis presented here, we set m = 10 milliseconds. 
This assumption also improves the computational prop- 
erties of our posterior simulations by preventing them 
from becoming "stuck" at values of the r parameter that 
are near the ends of the light curve. The end result of 
these restrictions is that our prior on r is uniform on the 
set [m, T — m].^ 

We can then use a Markov chain Monte Carlo (MCMC) 
technique (the Metropolis-Hastings algorithm) to draw 
from the posterior distribution of our parameters. From 
these draws, we can calculate quantities such as the pos- 
terior median and 68% posterior intervals for (hi — /12) 
and (fci — to characterize the deviation of our events 
from symmetry and our uncertainty about these prop- 
erties. For an excellent introduction to the use of the 
Metropolis-Hastings algorithm for statistical inference, 
please consult Chib & Grccnbcrg (1995). 

3. COMPUTATIONAL DETAILS OF POSTERIOR 
SIMULATIONS 

Although the model presented above is quite appeal- 
ing, the resulting posterior distribution for our param- 
eters of interest is not analytically tractable. Thus, to 
make inferences about these parameters, we simulate 
a series of samples from their joint posterior distribu- 
tion. We use a blocked Metropolis-Hastings algorithm 
for these simulations. This means that, instead of draw- 
ing new values for all seven of our parameters at once 
and accepting or rejecting them based on the Metropolis- 
Hastings rule, we vary only a subset of the parameters at 
time, holding all others fixed. The blocks used are: (1) r, 
(2) c and a, and (3) hi, /i2, ki and fc2. So, for each itera- 
tion of our algorithm, we first draw a new value of t from 
our proposal distribution and choose whether to accept 
it based on the Metropolis-Hastings acceptance probabil- 
ity. This acceptance probability is calculated holding the 
values of all other parameters (c, a, hi, h2, etc.) fixed. 
Next, we draw new values of c and a from our proposal 
distribution for these parameters and calculate an anal- 
ogous acceptance probability. Finally, we perform a sim- 
ilar procedure on the shape parameters {hi, h2, ki, and 
fc2). Blocking in this way allows us to take advantage of 
the structure in the model in constructing our simulation 
algorithm, improving its efficiency. For more information 
on blocked Metropolis-Hastings methods, please consult 
§6.2 of , Chib fc Grcc nbcrg (1995). 

3.1. Initialization 

We look to the data to obtain starting values for our 
MCMC algorithm. The first step is to obtain a smoothed 
version of our signal yt, which we will denote yt- This 
is generated by taking a centered moving average of the 
original signal yt- For the analyses presented here, a 
window ten milliseconds wide was used, although this 

^ However, for other analyses of time symmetry in events (such 
as the analysis of longer duration events that are not fully contained 
in a light curve), these assumptions could be relaxed. 



can be modified. The important point on the selection of 
the window is that it should be narrow enough to localize 
the event while still performing some noise reduction. If 
the window is too wide, our smoothed light curve will 
not capture the structure of the the event; as a result, it 
would not be useful in determining initial values for the 
magnitude (a) or in determing the location of the event 
in the light curve. However, if the window is too narrow, 
it will not smooth the light curve well. As a result, our 
initializations would be far noisier than necessary. 

Using yt, we can construct our initialization values and 
a proposal distribution for t. Focusing on initialization 
first, we set cq = median (yt) and ao = Cq — mint(yt). 
We use a robust measure of typical baseline rate (the 
median). The motivation for this initialization of c is 
that, if the event is brief relative to the duration of the 
light curve and our smoothing window is not too long, 
our smoothed light curve should be near the baseline 
count rate c for most times. In contrast, if we simply 
used the overall mean, it would be contaminated by the 
event of interest. Our initialization for a is relatively 
straightforward; we take our estimate of the baseline rate 
(co) and subtract the minimum observed count rate to 
obtain an estimate of the magnitude of our event. 

To initialize hi and /12, we take a similar approach. 
First, we assume that a symmetric initialization will be 
an acceptable starting point for our analysis. We then 

calculate ho ^ ^# {t : yt < "'°dian (^,)+min,(i..) | 

hi and /12 to hg.^ Based on the above discussion of co 
and ao, this approach to initializing hi and h2 appears 
quite natural. We are using our initial estimates of c and 
a to estimate the FWHM of our event and attributing 
equal parts of this time to hi and /12. Finally, we initial- 
ize ki and fc2 manually. Our posterior simulations were 
quite insensitive to the initializations of these parame- 
ters; values between 1 and 20 yielded similar results. 

3.2. Metropolis-Hastings steps 

Keeping in mind the prior discussed in |§2.2[ we can 
obtain a proposal distribution for t using the procedure 
outlined in Algorithm 1. 



Algorithm 1 Construction of proposal distribution for 

T 

Require: yt for i G {1, . . . , T}, mGN, p>0, e>0 

wt ^ (maxt (yt) - ytY 

wt ^ wt + t 

if t < m or {T — t) < m then 

Wt ^ 
end if 



Calculate points along cumulative distribution function (CDF) 
as Pt = T,s<t 

Calculate CDF F{t) via linear interpolation of {pt} 
Calculate quantile function q{p) via linear interpolation 
Calculate density function /(i) as derivative of F{t) for t ^ N 



At the conclusion of this procedure, we have the den- 
sity, cumulative density, and quantile functions for our 
proposal distribution for r. Through a series of numeri- 
cal experiments, we found that values of p between two 

® # is the cardinality operator, which gives the number of items 
in the given set 
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and five typically produced the best proposal distribu- 
tions; p was set to four in the subsequent analyses. The 
e term serves to ensure that all Wt are non-zero, avoiding 
an issue with the interpolation of our quantile function 
(where a wt of zero leads to a vertical segment). We use 
linear interpolation because it is the only spline method 
that produces identical functions when the axes are ex- 
changed. This is vital in our application because we are 
conducting seperate interpolations of {pt} to obtain our 
CDF and its inverse (the quantile function). If we had 
used any other type of polynomial interpolation, the in- 
verse of the calculated spline for the CDF would not be 
another polynomial spline. However, with a linear spline, 
the inverse of the interpolated function is simply another 
linear spline. Thus, we can avoid complex calculations 
and use standard interpolation routines to calculate our 
quantile function. Our method of generating a proposal 
for our location parameter is quite versatile, adapting 
naturally to the presence of multiple potential events in 
the light curve. We will discuss some potential applica- 



tions of it beyond this analysis in § 6 

The proposal generated by the above procedure allows 
us to use an independence chain Metropolis-Hastings 
method for r. With this class of methods, the proposal 
distribution for the parameter being simulated is inde- 
pendent of that parameter's current value in the MCMC 
simulation. If the proposal is well-constructed, this can 
greatly improve the convergence properties of the simula- 
tions; in the ideal case, we would make our proposal the 
same as our posterior to maximize efficiency. However, 
for such a procedure to be effective, one must have some 
prior knowledge of the parameters 's posterior distribu- 
tion. That is what the above procedure extracts from 
the light curve. 

Using the initial values and proposal distribution ob- 
tained above, we construct proposal distributions for the 
parameters of our model (excepting r) . All proposals are 
Gaussian and are thus characterized by their means and 
standard deviations: 



(6) 

■) (7) 

) (8) 



(c,a) : ^1 = (co,q;o), a = {\j — , y/a^) 
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With these proposals and the blocking scheme de- 
scribed previously, we have fully specified our MCMC 
simulation approach. We now demonstrate it's effective- 
ness with simulated data. 

4. RESULTS ON SIMULATED SYMMETRIC EVENTS 

To test our approach, we simulated a 500 time steps 
light curve with a symmetric event resembling those ob- 
served in the Scorpius X-1 data. For this simulation, 
we set c = 50, a = 40, hi = /i2 = 5, fci = ^2 = 2, 
and r = 300. The resulting light curve can be seen in 
[Figure 3] 

The r esults of the poste rior simulatio ns are summa- 
rized in [Figure 4[ [Figure 5[ and [Figure 6[ From the first 
of these ( [Figure 4p , we observe that the posterior dis- 
tributions loFToth the baseline count rate (c) and the 
magnitude of the event (a) are centered around their 



true values of 50 and 40, respectively. The posterior on c 
is significantly more concentrated than the posterior on 
a, as we would expect; we can use nearly every point in 
the light curve to infer the value of c, but we can only use 
the points during the event to infer the value of a. The 
posterior distributions of (7z4_^_/i2) and (fci — are con- 
siderably more diffuse, as [Figure 5[ demonstrates. These 
posteriors are, however, centered near the true value of 
0. The relationship between the posterior distributions 
of these two quantities can be seen from [Figure 6[ There 
appears to be a significant negative dependence between 
the two deviations. 

The properties of our estimates for symmetry-related 
parameters ((/ii — ft.2) and (^1 — ^2)) are further supported 
by a second set of simulations. Using the parameter val- 
ues given above, we simulated 1000 light curves and ran 
our algorithm on each one. This allows us to quantify 
some frequentist properties of our Bayesian estimators, 
such as the coverage of our posterior intervals.'' The re- 
sults of these simulations are summarized in [Figure T\ 
Both distributions appear centered around zero! Cov- 
erage properties are not perfectly calibrated (as this is 
a Bayesian method), but they appear reasonable. 938 
of our 1000 68% posterior intervals for {hi — ^,2) in- 
cluded zero, as did 816 of our 68% posterior intervals 
for {ki - k2). 

Overall, the results of this exercise with simulated 
events serve to increase our confidence in the validity 
of our method. We now apply our model and algorithm 
to the analysis of the Scorpius X-1 dip events. 

5. RESULTS ON SCORPIUS X-1 EVENTS 

Usi ng the above model, we analyzed all 107 events 
from 'Cha ng et all (|2007f ). For each event, we ran our 
MCMC algorithm for 400, 000 iterations, discarding the 
first 200,000 as burn-in^. Run in parallel on Harvard's 
Odyssey cluster, these runs required approximately 400 
seconds in total (approximately 43,000 CPU seconds). 
Typical acceptance rates for the individual Metropolis- 
Hastings steps were between 10% and 40%. The re- 
sults of t hese simulation s are summarized in [Figure §1 
[Figure 9l and [Figure 10[ 

From [Figure 8[ we see that the FWHM of a typical 
event from .Chang et ahl (|2007[ ) is between 2 and 4 mil- 
liseconds (with 68% posterior probability) . Events 2 and 
107 appears to be the only significant deviation from this 
pattern with FWHMs of approximately 7 and 13 millisec- 
onds, respectively. Visual inspection of these events (in 
Figures [H] and (T^l) reveals that both are unusual in ap- 
pearance. Event 2 is characterized by a more gradual 
change in intensity than most and appears to have two 
minimia. Event 107 appears to have two clear, significant 



ap- 
ow- 



lliiiiiillU. 

The posterior intervals for [hi — /i2) in [Figure 9 
pear to be relatively evenly scattered around zero. E 

^ The coverage of an interval estimator is defined as the percent- 
age of intervals that will cover the true value of the parameter of 
interest across repeated sampling. By checking the coverage of our 
key estimators, we can quantify their classical error rates in addi- 
tion t o their Bayesian properties. For further details, Wassermanl 
pool ) provides an excellent reference in §6.3.2 and §11.9. 

* In MCMC simulations, it is common to discard the earliest 
portion of the results as burn-in to ensure that we are only bas- 
ing our conclusions on draws from the Markov chain's stationary 
distribution (which is, by design, our posterior). 
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ever, although most of these intervals (90 of 107) cover 
zero, there is notable tendency towards negative values 
of {hi — h2)- This can be seen in the posterior median, 
with 81 of the 107 posterior medians below zero. Using 
a classical sign test of the hypothesis that each poste- 
rior median is equally likely to be positive or negative 
pixon fc MooJl94d) . we obtain a p- value of 9.4 x 10"*^, 
indicating that there is a significant deviation from sym- 
metry in the relative half-maxima for this collection of 
events. The sign test is conducted by first counting the 
number of positive posterior medians. We then calcu- 
lated the probability of observing a result at least as 
far from an even split of positive and negative medians 
as was observed, assuming that each median is equally 
likely to fall above or below zero. In this case, the p- 
value was calculated as 2 x P{X < 26\p = where X 
is the number of observed positive medians. Under the 
null hypothesis that positives and negatives are equally 
likely, X ~ Bin(107, i), so the preceding probability can 
be easily calculated using standard techniques for the bi- 
nomial distribution. The pattern of negative deviations 
in (fci — can also be seen in [Figure 13} the devia- 
tio n of this distribution from the symmetric case (seen 
in [Figure 7"| is cle ar. Event s 2 and 107 appear some- 



what anomalous in 
unusual structure o: 



Figure 9| this is consistent with the 
' these events discussed above. 
Finally, the posterior medians for (fci — ^2) are con- 
centrated below zero (90 of 10 7 are nega tive). Inspec- 
tion of the posterior intervals in [Figure 10[ indicates that 
very little posterior support is typically found above zero 
(only 55 of the 107 intervals included zero). A classical 
sign test on the posterior medians (as above) yields a 
p-value of 3.5 x 10~^^. This indicates that the observed 
pattern of posterior medians is quite unlikely to occur if 
the true distribution of posterior medians is evenly dis- 
tributed above and below zero. This is consistent with 
the findings of ILiu et al.l (j2008l ). as a negative value of 
(fci — ^2) indicates that the count rate decreased more 
quickly than it increased. The pattern of negative de- 
viations in (fci — ^2) can also be seen in [Figure 14[ the 
deviation of this distribution from the symmetric case 
(seen in [Figure 7[ ) is again apparent. 

The concentration of {hi — ^,2) and (fci — across 
events indicates that the observed events are character- 
ized by a relatively rapid decrease in the count rate fol- 
lowed by a more gradual return to the baseline level; that 
is, the typical event profile is right-skewed. This pattern 
can be seen in |Figurc 15 which shows four representative 
events from lChang et af] ([20071 ). 

It is interesting to note that {hi — ^,2) and (fci — fc2) 
are typically negatively correlated (when examining their 
joint posterior distribution) for approximately symmet- 
ric events. An example of th i s com es from our analysis 
of event 90 from lChang et al.l (|2007D . The joint posterior 
distribution of { hi — /12) a nd (fci — fc2) can be seen (as 
a scatterplot) in [Figure 16[ Two points are immediately 
apparent from this hgure. First, the deviation of the 
distribution from the origin indicates that this event is 
most likely asymmetric. Second, {hi — /12) and (fci — fc2) 
are strongly negatively correlated (r — —0.618). The in- 
tuition behind this outcome is clear. We are allowing 
for two possible deviations from symmetry in our model: 
deviations based on the half maxima {hi — /12) and devia- 



tions based on the rate of change (fci — fc2) in our source's 
intensity. For very short events, we cannot be certain to 
which of these deviations we should attribute the asym- 
metry of our event, as we can, to some extent, trade-off 
between them. Our posterior distribution simply reflects 
this trade-off. 

It is also important to note that, when analyzed us- 
ing our method, the events identified as potentially non- 
instrumental in origin by [Chang et al.l do not appear 
to deviate significantly from the remaining events in 
FWHM, half maxima deviation, or curvatu re deviation. 
This supports the conclusion of ILiu et all |2P08) that 
the events detected thus far (including those identified 
as possibly non-instrumental) are most likely the result 
of dead-time effects and other instrumental contamina- 
tion, not TNOs occulting Scorpius X-1. 

At the distance of the Kuipcr Belt (^ 40AU) the size 
of the objects given the width of the events is ~ 50 m, 
and that is close to the Fresnel scale, which is -^Z A a/2, 
where A is the wavelength and d is the distance. We use 
A = 0.3 nm (4 keV) because most of the RXTE/PCA- 
detected photons from Sco X-1 are at this energy. As a 
result these occul tation events a r e diff ractio n dominated 
pheno mena (see iRoques et al.l (|1987[ ) and iNihei et "all 
(|2007l )) and we expect that diffraction effects would pro- 
duce small count rate increases on either side of the dips. 

Also TNOs of t his size are not exp ected to be spher- 
ically symmetric ([Rogues et al.l ll987f) and therefore the 
resulting occultation events in the light curves are not ex- 
pected to be time symmetric. However for TNOs of the 
size of the Fresnel scale or smaller the shape of the ob- 
ject does not translate into asymmetry of the time series. 
Roques et al 1986 has studied light curves from different 
shaped TNOs and demonstrated that when the effective 
size of the TNOs is smaller than the Fresnel scale the re- 
sulting events in the light curve do not differ from the 
light curves signatures of regularly spherically shaped 
TNOs. The sizes of the 107 objects studied here are 
within the range of the Fresnel scale so some asymmetry 
on the light curves should be manifested. Figures [TTl and 
118! show the measures of asymmetry as a fu nction of the 
size of the object for 57 of the 58 events from lChang et all 
(|200l) . As can be seen from the figures there is no clear 
relation between the estimated size of an object and the 
level of asymmetry for its corresponding event. This pro- 
vides further evidence that the observed dips in intensity 
do not correspond to TNO transit events, as 1 Jones et all 
([2008D suggest. 

6. CONCLUSIONS 

Based on our analysis of the suspected Scorpius X-1 
occultation events, we conclude that there is significant 
evidence for temporal asymmetry in these events. This 
findin g is consistent with the suggestion of iJones et al.l 
(|2008l ) that the observed dips in brightness are due to 
dead ti me effects. We als o conclude that the events iden- 
tified in ILiu et al.l (|2008f ) as possibly noninstrumental in 
origin are most likely not occultations. 

We also note that the approach presented herein can 
be applied (with minor modifications) to the study of 
events in nearly any light curve. Such studies need not 
be limited to examining the symmetry of events, as the 
parameterization of our bilogistic event profile and flex- 
ible Bayesian approach allow inference to be made on 
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quantities such as full width at half maxima. For ex- 
ample, our approach could be extended to analysis of 
exoplanetary transits, supernovae, quasars or any other 
transient event. 

Furthermore, the method we have presented for con- 
structing a proposal distribution for the center of each 
event can be applied quite generally. Even when us- 
ing a physically grounded model with parameters such 
as body size, angle of incidence, etc., one must account 
for uncertainties about the location of the event in the 
light curve when making inferences about the parame- 
ters of interest. This can be quite challenging using clas- 
sical methods, and even MCMC methods do not offer an 



panacea. Our approach to constructing an independence 
chain proposal distribution for the location parameter 
in such models makes them quite a bit more tractable, 
especially in the case of multiple events. 
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Fig. 1. — : Bilogistic event profile for c = 50, a = 40, fci = ^2 = 1, and r = 50. Varying hi and /i2 from 1 to 16. 
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Fig. 2. — : Bilogistic event profile for c = 50, a = 40, /ii = /12 = 8, and r = 50. Varying ki and ^2 from 1 to 8. 
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Fig. 4. — : Histogram of posterior distributions of c and a. for simulated event 
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Fig. 5. — : Histogram of key symmetry-related parameters for simulated event 
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Fig. 7. — : Distributions of posterior medians of (A;i — and {hi — /12) for 1000 simulated events 
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Fig. 8. — : 68% posterior Intervals & medians for FWHM {hi + ^,2) for IChang et al.l ()2007n occultation events. 
Second plot is zoomed. 
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Fig. 9. — : 68% posterior Intervals & medians for {hi — /12) for lChang et alJ (|2007f ) occultation events. Second 
plot is zoomed. 
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Fig. 11. — : Light curve for event 2 from lChang et all (|2007f ). 
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Fig. 13. — : Histogram of posterior medians for {hi — h2) for lChang et al.l (|2007f ) occultation events. 
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Fig. 14. — : Histogram of posterior medians for {ki — k2) for lChang et al.l ()2007f ) occultation events. 
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Fig. 16.—: (fci - fe) vs. {hi - ^12) for event 90 from lChang et all (|2007f) . 
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Fig. 17. — : \hi — /i2| vs. estimated diameter for 57 of the 58 events from I Chang et all (|2006[ ) (58th event omitted due 
to anomalous nature). 
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Fig. 18. — : \ki — vs. estimated diameter for 57 of the 58 events from lChang et all (|2006f ) (58th event omitted due 
to anomalous nature). 




